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Abstract 

Relations connecting various zero momentum correlators of interpolating 

fields for pseudoscalar and scalar channel, containing one heavy and one light 

quark field, are derived from the Euclidean space formulation of the QCD 

functional integral. These relations may serve as constraints on the phe- 

nomenological models or approaches motivated from QCD, and suggest a 

method to extract the chiral quark condensates. It is also found that the 

correlator for pseudoscalar channel differs from that for scalar channel even 

in the large heavy quark mass limit. 
PACS numbers: 11.15.Tk, 12.38.Aw 
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Typeset using REVT^ 



Heavy-light systems made of one heavy and one hght quark have been attracting atten- 
tions, motivated by the observation of heavy quark symmetry in the infinite heavy quark 
mass hmit fl]]. The basis of this symmetry is a separation between the scale of the heavy 
quark mass and all other dynamical scales in the problem. 

The properties of heavy-light systems have been investigated using various phenomeno- 
logical approaches such as QCD sum-rule method and numerical simulations on a lattice 
IBHlOi, based on the analytic properties of two-point correlators of interpolating fields, car- 
rying the quantum numbers of the system under study, and on asymptotic freedom. In these 
approaches, hadronic spectral properties {e.g., masses, coupling constants, etc.) are related 
via dispersion relations to the correlators evaluated in terms of quark and gluon degrees of 
freedom. In practical applications, one needs to parametrize the spectral functions and to 
evaluate the correlators approximately. The success of such approaches depends on correct 
understanding of the qualitative features of the spectral functions and accurate evaluation 
of the correlators from QCD. 

In this paper, we show, within the Euclidean formulation of the QCD functional integral, 
that the two-point zero momentum correlators of interpolating fields for both pseudoscalar 
and scalar channel containing one heavy and one light flavor can be expressed as combi- 
nations of quark condensates. We then derive relations which connect various correlators. 
These relations can be regarded as constraints on phenomenological descriptions of the cor- 
relators. We find that the correlators for pseudoscalar and scalar channel are different even 
in the large heavy quark mass limit, suggesting differences between positive- and negative- 
parity states. 

To evaluate various two-point correlators from QCD, it is necessary to know the proper- 
ties of the QCD vacuum. This is, however, difficult due to the nonperturbative features of 
QCD at large distances, which lead to complicated properties of the QCD vacuum such as the 
formation of the gluon condensate (0|(as/7r)G'^|0) and the quark condensate (0|gg|0). The 
latter is of particular importance as it signals spontaneous breaking of the chiral symmetry 
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which the QCD Lagrangian enjoys in the massless hmit. 



In the pioneering work of Banks and Casher , it was noted that the fermion condensate 
is closely related to the mean level density of the fermion eigenvalues (0|gg|0) ~ —np(\ = 0), 
where p{X)dX denotes the mean number of eigenvalues of the Euclidean Dirac operator, p, 
contained in the interval dX, per unit volume. Thus the formation of a quark condensate is 
connected with the occurrence of the unusual behavior of the low eigenvalues of the Dirac 
operator, which is dramatically different from that for the free Dirac operator. 

The Bank-Casher formula (0|gg|0) ^ — 7rp(0) can be easily derived in the Euclidean 
formulation of the QCD functional integral |T2|-p!^. In a given background gluon field 



configuration, the fermion Green's function can be written as |12 



Giix, y) = {qf[x)qf[y))A = 2^ tt— , (1) 

„ rrif — tAn 

where / is flavor index (a sum over / is not implied), and ipni^) and A„ are normalized 
eigenfunctions and eigenvalues of the Euclidean Dirac operator: 

PM^) = XnM^) ■ (2) 

Note that except for zero modes (A^ = 0), the eigenfunctions occur in pairs, {ipn, TsV'n}, 
with opposite eigenvalues, {A„, — A„} [0,0. Consequently, one obtains 

iXA.T.{e(,...)}4E;|:%. (3) 

where V is the volume of the four- dimensional box introduced to provide the infrared regu- 
larization of the theory and make the spectrum discrete; at the end, one should take the limit 
V oo. The zero-mode contributions, vanishing in the infinite volume limit, have been 



droped in Eq. @ (see Ref. [jT2[ for details). Averaging over all gluon field configurations, 
one gets for large volumes 

{Olq^qjlO) = ^2mf f dX^^^^ , (4) 
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where p{X) is the mean level density. Here the average means the functional integral over 
gluon configuration with weighting function, e""^"^^ mutiplyed by the fermion determinant 
functional, where Syu is the Euclidean action for pure Yang-Mills theory. Note that the 
integral in Eq. (^) is divergent due to the perturbative ultraviolet contributions. One can, for 
example, regularize this divergence by subtracting the perturbative contributions from p(A) 
or simulate the functional integral in a lattice regularized version. In the discussions to follow, 
we always assume that the integrals over the mean eigenvalues are properly regularized. Our 
results are independent of how the associated divergences are regularized. Assuming that 
the quark mass is much less than the characteristic hadronic scale on which p(A) changes 
substantially, one immediately reproduces the Bank-Casher formula. 

To begin, we consider the correlators of pseudoscalar interpolating fields containing one 
heavy and one light quark field as given by 

n^f (?) = ^ / d'xe^''-{0\TfJ{x){fjymO) , (5) 

where the colorless interpolating field j^f is defined as 

jhh^) = i^i + mh)l{x)i-f5h{x) , (6) 

where l{x) denotes light quark field {u, d, s) with mass mi and h{x) is for heavy quark 
field (c, b, t) with mass nih. The interpolating field j^f has the same quantum numbers 
as a pseudoscalar heavy-light meson {D, B, etc.), and is equal to the divergence of the 
axial- vector current, A'/^i{x) = I(x)7^75/i(x). Since is partially conserved, U^f{q) is a 
renormalization-group invariant quantity. In addition, Lorentz covariance, parity and time- 
reversal imply that Il^f is a scalar. 

Calculation of li^f{q) from QCD is in general difficult. However, at zero momentum, 
the correlator H^f (0) can be easily evaluated within the Euclidean formulation of the QCD 
functional integral. For a fixed gluon configuration, only the connected part contributes, 
and one gets 
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= {mi + m^f^ ^ ATr [qI^x, y)i^Q\{y, x)^^ , (7) 
where Q^p^ and Q\ have been defined in Eq. (|T]). A simple alegebraic calculation then yields 



nPf(O) |a= K + m,)- ^ 



An>0 



2mi ^ 2m/i 



(8) 



Averaging over the gluon background and using Eq. (H), we obtain 



(0) = + mn) (0|I/|0) + (0|M|0) 



(9) 



Thus, the zero momentum correlator IlJ^f (0) is related to the quark masses and quark con- 
densates. The analytic properties of IlJ^f (0) allow one to write a dispersion representation, 
where the spectral function contains the information about all the physical states coupling 
to the interpolating field j^f . If IlJ^'f (0) were saturated by the lowest physical resonance, 
one would expect to obtain a relation from Eq. (^), similar to the Gell-Mann-Oakes-Renner 
relation, which relates the properties of the lowest resonance to the quark masses and quark 
condensates. In nature, however, the heavy quark masses are much greater than the QCD 
scale parameter, and the correlator IlJ^f (0) is not saturated by the lowest resonance. 

Nevertheless, some useful relations can be immediately derived from Eq. (P). For two 
different light flavors and a flxed heavy flavor, we have 



n^f (0) - n;:,^(o) = (o|/'/'|o) - (o|I/|o) 



rPS 



(10) 



where I and /' are different light flavors and h can be any heavy flavor. For two different 
heavy flavors and a flxed light flavor, we get 



nCf (0) - — nP^(o) = [(o|/i'/.'|o) - (o|/./i|o) 



(11) 



where h and h! are two different heavy flaovrs and / can be any light flavor. Finally, for two 
different heavy flavors and two different light flavors, we flnd 
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n]:f (0) - ul!{0) - ^ul^io) + ^ulUo) = , (12) 

nih' rrih' 

where h and h' are different heavy ffaovrs and / and /' are different hght ffavors. In Eqs. (p!0|- 
|I2|), we have assumed that rrihy ^ nii^ii and used nih^h' ± — 

Note that in deriving Eqs. (p|-p!^) we have not made any assumptions, except for nihy 3> 
mi J/. So, they are all model independent results. (If one retains rrih^h' explicitly in the 

derivation, the resulting relations will be exact). If one adopts dispersion representations 
for the correlators, these results will constrain the spectral functions. Phenomenological 
parametrizations of the spectral functions, often used in applications such as the lattice 
QCD data or QCD sum-rule calculations, must satisfy these constraints. 

One notices from Eq. (|ToD that 

(o|Sd|o) - (o|m|o) = [nL^(o) - nL^(o)] , (i3) 
{o\ss\o) - (o|m|o) = [nL^(o) - nL^(o)] . (i4) 

Therefore, one may determine the SU(3) breaking in light quark condensates by evaluating 
the correlators Ul^iO), 0^1(0), and H^f (0). 

Following the same pattern as used in deriving Eq. (H), one can easily obtain 

nPS(o) = -Am,{0\hh\0) , (15) 

where 11^^ is the correlator of interpolating field j^hi^) — '^'mhh{x)i'~^^h[x) . Combinig this 
result with Eq. (||), we have 

(o|^«|o) = — [nPS(o)-4nL^(o)] . (i6) 

So, the calculation of HIH^) and H^f (0) will give rise to a determination of the quark con- 
densate (0|mm|0). In the lattice calculations, one usually determine the quark condensates by 
extracting p(A = 0) directly [|T5l -[T7[|. It will be interesting to extract the quark condensates 
by calculating various correlators discussed here on a lattice, and compare the results with 
those obtained from the direct calculation of the mean level density. 
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We now turn to the scalar channel. The corresponding correlator is defined by 



UUq) ^ z J rfW^-(0|Tjl,(x)(jl,)t(0)|0) , (17) 

where jhi{x) = {mi + mh)l{x)h{x) is the scalar interpolating field. Again, we are interested 
in zero momentum correlator 11^; (0), which, in a fixed background gluon configuration, can 
be written as 



(18) 



Taking into account the pairing properties of the eigenfunctions and completing the integral, 

we get 



/ N , / n9 1 2A^ — 2mimh 

n?.,(o)U^(™, + ™OV£KTI|)KTA;) 



A„>0 



{mi + TTihY 1 
mh-mi V 



E 

A„>0 



2mh 



2m; 



(19) 



Averaging over all gluon configurations and using Eq. (^), one obtains 



(0|I/|0) - {0\hh\0) 



(20) 



rrih - mi 

Similarly, one can derive relations which connect different correlators. Applying Eq. (^OD to 
two distinct light fiavors, one gets 



n^,(o)-n^,(o) = m, (o|I/|o) - (o|/'/'|o) 



For a fixed light fiavor and two different heavy flavors, one obtains 



— n^.(O) - UUO) = m, [{0\hh\0) - {0\h'h'\0) 
my L 



With two different light flavors and two different heavy flavors, one then flnds 

n^.(o) - nL(o) - ^nlm + ^nl,{o) = o . 

rrih' my 



(21) 



(22) 



(23) 



Here we have neglected light quark masses (i.e., mh^y ± m/ ~ mh^y). These equations are 
model independent, and can be regarded as constraints on phenomenological descriptions of 



the correlators for scalar channel. One may also use Eqs. (pOD and (|2T|) to extract the light 
quark condensates. 

Combining Eq. (P) and Eq. (pO]), we find the following relations 

Ulfi0)-Ul{0) = -2m,{0\m, (24) 

nPf(o) + n^,(o) = -2m,(o|Ho), (25) 

which connect the correlator for pseudoscalar channel to that for scalar channel. One notices 
that these two equations can be used to determine both light and heavy quark condensates. 
These relations also show that even in the limit rrih oo, the pseudoscalar channel differs 
from the scalar channel. It is worth emphasizing that this result does not contradict the 
heavy quark symmetry. In the infinite heavy quark mass limit, heavy quark symmetry 
implies that the spin and mass of a heavy quark decouple from the hadronic dynamics; 
consequently, there are degeneracy in spin, but not in parity. 

From Eqs. (IT^); (O); (0); (P^), one may also deduce the following relations 



n^f (0) - nP^(o) + n^,(o) - n^,(o) = o , (26) 
mo) - nL(o) - ^n^,^(o) + ^n^,(o) = o . (27) 

In conclusion, we have derived exact relations which connect various zero momentum 
correlators of interpolating fields made of one heavy and one light quark field for pseudoscalar 
and scalar channel, within Euclidean formulation of QCD functional integral. These relations 
may serve as constrains on phenomenological models for the correlators such as those used 
in QCD sum-rule calculations. We also noted that in principle one can use these relations 
to extract quark condensates by calculating the correlators on a lattice. One can extend the 
present techniques to study the vector and axial vector channel. 
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